Two-photon cooling of a nonlinear quantum 

oscillator 



Mihai A. Macovei ^ , Xiao-Tao Xie 

Max-Planck-Institut fiir Kernphysik, Saupfercheckweg 1, D-69117 Heidelberg, 

Germany 



Abstract 

The cooling effects of a nonlinear quantum oscillator via its interaction with an 
artificial atom (qubit) are investigated. The quantum dissipations through the en- 
vironmental reservoir of the nonlinear oscillator are included, taking into account 
the nonlinearity of the qubit-oscillator interaction. For appropriate bath temper- 
atures and the resonator's quality factors, we demonstrate effective cooling below 
the thermal background. As the photon coherence functions behave differently for 
even and odd photon number states, we describe a mechanism distinguishing those 
states. The analytical formalism developed is general and can be applied to a wide 
range of systems. 
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1 Introduction 



Simple models describing the main properties of various phenomena in physics 
are always of particular importance. The Jaynes-Cummings model, for in- 
stance, provides a simple description of the interaction of matter with an 
electromagnetic field pQ. It consists of a two-state particle interacting with 
a single quantized mode, applicable in general to cavity quantum electrody- 
namics. Now, there is an increased interest to apply such a simple model to 
more complex systems like superconducting electrical quantum circuits. This 
allows us to investigate them in an analogous way as a two-level atom inter- 
acting with a quantized electromagnetic cavity mode [2|. As an advantage, for 
instance, the strong coupling limit and enhanced lifetimes can be achieved in 



^ On leave from Institute of Applied Physics, Chi§inau, Moldova 

Email address: mihai.macovei@mpi-hd.mpg.de (Mihai A. Macovei). 



Preprint submitted to Elsevier 



22 March 2010 



superconducting devices [3]. Remarkably, a Josephson qubit was designed that 
entangles qubit information P] . Entanglement between a superconducting flux 
qubit and a superconducting quantum interference device was demonstrated 
in [S], while a procedure to directly measure the state of an electromagnetic 
field inside a superconducting transmission line, coupled to a Cooper-pair 
box, was proposed in [6], making them attractive to quantum computation 
processing. Of particular interest are the studies regarding the dephasing of 
the superconduction qubit induced by the photon noise [7|. Other systems re- 
fer to coupling of the superconducting qubit to a solid-state nanomechanical 
resonator. Various interesting effects related to the nanomechanical resonator 
state were demonstrated. In particular, squeezing of the nanomechanical res- 
onator state occurs when coupling it to a Josephson quantum circuit [8]. The 
fidelity of a state transfer from the Josephson junction to a nanomechani- 
cal resonator was investigated in [9]. Entanglement from a nanomechanical 
resonator weakly coupled to a single Cooper-pair box [10], continuous mea- 
surement of the energy eigenstates of a nanomechanical resonator without a 
nondemolition probe [11] or signatures for a classical to quantum transition 
of a driven nonlinear nanomechanical resonator [12] were already discussed. 

Via engineering superconducting elements as artificial atoms and coupling 
them to a photon field of a resonator or to vibrational states of a nanome- 
chanical resonator one can demonstrate other interesting phenomena such 
as single artificial atom lasing or cooling. In particular, schemes for ground- 
state cooling of mechanical resonators were proposed in [13\. A flux qubit 
was experimentally cooled [2] by using techniques somewhat related to the 
well-known optical sideband cooling methods (see e.g. Ref. [15] and references 
therein). Continuous monitoring of Rabi oscillations in a Josephson flux qubit 
was reported in [16J while lasing effects of a Josephson-junction charge qubit, 
embedded in a superconducting resonator, was experimentally demonstrated 
in [T7] . Single-qubit lasing and cooling at the Rabi frequency was proposed in 
|18] . while a mechanism of simultaneously cooling of an artificial atom and its 
neighboring quantum system was analyzed in [19]. In some of these systems 
the nonlinear qubit-oscillator interaction was considered, i.e. the case when 
the qubit exchanges simultaneously more than one photon with the resonator 
mode. However, the dissipations of the nonlinear system due to interaction of 
the quantized mode with the environmental reservoir are more complex and 
require a further treatment. 

Thus here we report additional results regarding the nonlinear matter-light 
interactions which are general and applicable to a wide range of systems. 
To this end, we investigate the properties of a quantum oscillator coupled 
nonlinearly with a driven qubit through two-photon effects and damped via 
the thermal environmental reservoir, and focus on cooling phenomena of the 
oscillator's degrees of freedom. Due to a high degree of correlations between the 
particles generated in such a two-photon process, we consider the nonlinear 
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damping of the generated photons. This allows us to describe the system 
by using the properties of su(l,l) algebra. In the steady-state we obtain a 
mean photon number well below unity and thermal limit. Consequently, an 
effective cooling mechanism via nonlinear processes is discussed. Further, we 
propose a scheme which is able to distinguish between even and odd photon 
number states corresponding to su(l,l) algebra via measuring the second- 
order photon coherence function (and/or higher-order photon correlations). 
In addition, photon statistics may show quantum features, i.e. an important 
step towards single-photon sources. 

The paper is organized as follows. In Sec. 2 we introduce the system of interest 
and derive the corresponding master equation. The next section 3 analyzes the 
results. Finally the summary is given in Sec. 4. 



2 The model 

Particularly, we consider a Josephson flux qubit coupled inductively to a slow 
LC oscillator. The frequency of the oscillator is much lower than the qubit's 
tunnel splitting, i.e. Uc ^ A. The qubit is driven with Rabi frequencies near 
resonance with the oscillator frequency that affect the oscillator, increasing 
its oscillation amplitude. Near the symmetry point (i.e. the energy bias e 
between the flux states is negligibly small) and after transformation to the 
qubit's eigenbasis, the Hamiltonian describing the systems is: 

H = Aaz/2 + Q cos {ujt)ax + Uca+a — gax{a + a+), (1) 

where the first term describes the qubit while the second one considers its 
driving by an applied AC magnetic flux with amplitude Q and frequency u. 
The last two terms describe the oscillator with frequency Uc = 1/ V LC as well 
as the qubit-oscillator interaction, respectively. Here g ~ MIpIco, where M is 
the mutual inductance, Ip the magnitude of the persistent current in the qubit, 

and Jco = \JuJc/2L the amplitude of the vacuum fluctuations of the current in 
the LC oscillator. a+ and a are the creation and annihilation operators corre- 
sponding to the oscillator degrees of freedom, while ai {i G {x, y, z}) are the 
Pauli matrices operating in the dressed flux basis of the qubit subsystem. As 
A ^ Uc, the transverse coupling in the Hamiltonian ([1]) is transformed into a 
second-order longitudinal coupling by employing a Schrieffer- Wolff type trans- 
formation, i.e. Us = exp (iS) with 5* = {g/A){a + a^)ay [TS1I2U] . By further 
using the rotating wave approximation with respect to u and diagonalizing the 
qubit term as well as applying the secular approximation, i.e. omitting terms 
oscillating with the generalized Rabi frequency, one arrives at the following 
Hamiltonian describing the nonlinear interaction between the qubit and the 
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oscillator: 



- go{aa+ + a+a)aj4:. (2) 

Here g2=2g'^ sin26/A gives the nonlinear qubit-oscillator coupling strength 
while go=4:g'^ cos 29 / A accounts for a frequency shift of the qubit's frequency. 
Further cot2^^ = 6ijj/Q, where 6co = A — u and where = ^ {6uy + Q"^ 
stands for the generalized Rabi frequency. The Hamiltonian ([2]) involves two- 
quantum processes, i.e. two-particle exchanges between the qubit and the non- 
linear oscillator, which means that the quanta are created and annihilated 
simultaneously in pairs. The particles generated via such a quadratic process 
are known to be highly correlated, i.e. a single photon pair behaves like a 
quasiparticle [2T| . 

The spontaneous emission damping of the qubit in this picture is given as 
Psp = -l^'\cT^,cT^p]- E 7^"^H^"S^°V]+H.c., (3) 

ai^a2e{+,-} 

where = rocos^^/2, 7^-) = rosin^^/2 and 7^°) = Tosm^29/8. 

The damping of the quantized oscillator mode depends on the environmental 
reservoir. In order to have a two-photon damping of the nonlinear oscillator 
we consider that the quantum oscillator couples with the environmental bath 
via the following Hamiltonian 



Hf = hvh% + 2hx{h^p- + P+h). (4) 

Here the operators &^(&) belong to the broadband reservoir of carrier frequency 
V and represent the photon generation (annihilation) operator for the bath. 
Such a reservoir can be obtained by assuming that the LC oscillator couples 
additionally with another circuit the frequency of which v is equal or close to 
2ujc- Eliminating the bath operators in the Born-Markov approximation one 
can arrive at the master equation describing the damping of the nonlinear 
oscillator. For further convenience we introduce the field operators 



(3+ = a\l2, (3- = aV2 and (3, = (a+a + l/2)/2 

which obey the commutation relations for su(l,l) algebra, i.e. [(3'^, (3~] = —2f3z 
and [(3z, (3"^] = ±/3^. These operators act on the corresponding bases states of 
the su(l,l) algebra in the following way: 
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/3+ 1 m) = \J{m+ l){m + 2j) | j, m + 1) , 



/3 \j,m) = ^Jm{m + 2j - l)\j,m- 1), 
/^^Ij,"^) = ("^ + j)lJ,"^)• 



(5) 



Here m G {0, 1, 2, ■ ■ ■ , oo}, while for a single mode field, as considered in our 
approach, the allowed value of the Bargmann index (i.e., j) is 1/4 (3/4) for 
an even (odd) photon number. The correspondence between the number state 
of the single mode field \n) and the su(l,l) basis states \j,m) is \n) <H- \j,m) 



We have derived the master equation corresponding to the damping of the 
nonlinear oscillator via two-photon processes which can be written as follows 
(see Appendix): 



with Hq = 2xn/3z — describing an additional shift of the oscillator mode 

frequency proportional to and the Lamb shift proportional to X; respec- 
tively, induced by the thermostat via an effective coupling constant x- Here k is 
the two-photon damping rate of the quantized mode while n is the mean ther- 
mal photon number at frequency 2uc- In fact, for n = one obtains the well- 
known nonlinear damping of a quantized cavity mode via two-quantum pro- 
cesses used in Cavity Quantum Electrodynamics (see for instance [23|24f25j ). 
Finally, the master equation characterizing our model reads as follows: 



where H = 2{u,+xn)Pz-xP^p- + {nR-2go(3,)aj2+g2{a-f3+ + p-a+), while 
Apsp and Apf are given by Eq. ([3]) and Eq. (Q, respectively. Note here that 
the form of the master equation ([7]) would be similar to the corresponding one 
describing a wide range of problems involving two-quantum processes as for 
instance the quantum dynamics of a single two-state particle or a collection 
of two-state particles possessing dipole-forbidden transitions, pumped with an 
intense laser field in two-photon resonance and damped at resonance via two- 
photon effects by an optical cavity containing a two-photon absorber. Similar 
nonlinear damping, as in Eq. dHj), can be applied to a cavity mode crossed by an 
excited flux of such dipole-forbidden emitters. The laser/maser phenomenon 
via two photons can be developed here as well. Other applications refer to 
quantum effects in the present scheme as, for example, the first- and second- 
order squeezing of the oscillator's quantum fluctuations. These studies will be 
presented elsewhere. 



for n = 2{m + j) - 1/2 p3]. 



Pf = -z[Ho,p] - K{l + n){[P+,rp] + [p/3+,r]} 




(6) 



p = -i[H, p] - Ap,p - Apf, 



(7) 
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3 Results and discussions 



A general analytical solution of Eq. ([7]) is not evident. However, one can obtain 
its solution for different regimes of interest, namely in the bad or good cavity 
limit. Therefore, we proceed by investigating the properties of Eq. ([7]) when 
the qubit's quantum dynamics is faster than that of the nonlinear quantum 
oscillator, i.e. in the good cavity limit. Below the photon saturation num- 
ber no = (T\\T±/2g^y^^, with = 47^°) + T\\ and T\\ = + 7^-), one 
can integrate the qubit's degrees of freedom to arrive at a master equation 
characterizing the quantum oscillator only: 



p - txW+f3-,p] = - + n) + r_) /3-p] 
-(Kn + r+)[/3-,/3+p] + H.c. (8) 

Here r± = gl{l ± (a,)o)/(2rx), with (a,)o = (7^-^ - 7^+^)/r|| being the qubit 
inversion in the absence of the resonator mode. The two-photon resonance was 
assumed, i.e. Qr — 2gQ{f3z) = 2{uc + xn), as well as the relation: k(1 + n) <^ 
92 < To- 

The steady-state solution for the diagonal elements of Eq. ([8]) is 



p, = Z ^exp[-a/3^], (9) 

where Z is determined by the requirement Tr(ps) = 1 and a = In 77, with 
T] = (k,(1 + n) + Tj^/(^Kn + r_|_^. The expectation values of the operators 
needed for evaluating the properties of the nonlinear oscillator are obtained 
from Eq. ([5]) and Eq. (Q. In particular, the nonlinear oscillator mean photon 
number, i.e. (n) = 2{Pz) — 1/2, and its second- and fourth-order correlations 
can be determined from the following expressions: 



(/3.)=J+ ^ 



1] — 1 
2(l+j(r/-l)) 
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(/^-'rv "'^^'";;!'-^"^^"» +^. (10) 

It can be observed here that when rj approaches unity, the result is a substantial 
increase in the photon number and photon correlations. This will lead to lasing 
instability phenomena so that Eq. (|8]) and its solution are not valid anymore. 
However, rj ^ 1 corresponds to the cooling of the nonlinear oscillator where 
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Fig. 1. (color online) The mean photon number of the nonlinear oscillator (n) as 
a function of detuning 5u}. The solid blue curve stands for n=4, the long-dashed 
line for n=2, while the short-dashed one corresponds to n=l. The solid green 
curve shows the saturation photon number uq. Here, a;c/27r=27.5MHz, K/27r=2kHz, 



A/27r=3GHz, g/ln =18MHz, ro/27r=0.5MHz, n = Jnj^ - (5a;)2, and j = 1/4. 



the application of solution ([9]) is justified below the photon saturation number 
no, because the mean photon number as well as second- and fourth-order 
photon correlations tend to lower values in this case. Note that the control 
parameter rj can be modified by adjusting the qubit's parameters as well as 
the detuning of the external driving field. 

Fig. ([T]) depicts the mean photon number in the nonlinear oscillator mode, i.e. 
(n) when j = 1/4, as a function of various parameters governing steady-state 
behaviors. As can be observed here, lower photon numbers can be achieved via 
a suitable choice of the parameters involved and below the thermal limit. By 
increasing the coupling coefficient g such that g2 <^ ^r, the cooling efficiency 
can be further improved. Evidently, the qubit is more in its ground dressed- 
state, i.e. {crz)o < {6u > 0), when the cooling occurs. On the other hand, 
inversion of the qubit population can be created via modifying the detuning 
6u, that is for 6uj < 0. Thus, the cooling of the nonlinear oscillator occurs 
when controlling the qubit's population quantum dynamics. Although we get 
lower photon numbers for the nonlinear oscillator mode, it will be not easy, in 
general, to achieve (n) ~ 0. Due to approximations used in our approach, we 
cannot increase the coupling g ((72) further since we have performed the rotat- 
ing wave approximation in the Hamiltonian ([2]). The counter-rotating terms 
have to be taken into account when proceeding to larger g {g2)- Neither in this 
case can the degrees of freedom related to the qubit's quantum dynamics be 
adiabatically eliminated because (72 ~ Tq. Other limiting factors may appear 
due to fluctuations of external parameters. However, improving the oscillator 
quality factor one can achieve better cooling in general. 

We focus further on the properties of photon coherences. The second-order 
coherence function, i.e. ^''^^•'(O) = 4(/9"'"/3~)/(n)^, and the fourth-order one, i.e. 



^(4)(o) = (/3+'/3-')/(/3+/3-)2, evaluated by using Eq. ^ and Eq. © 
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Fig. 2. (color online) The second-order photon coherence function (^(^^(O) as a func- 
tion of detuning duo and for j = 1/4 (solid line) , and j = 3/4 (long-dashed curve) , 
respectively. The other parameters are the same as in Fig. ([1]) with n = 2. 

and represented as follows: 



(2).n^_ 32(l + (r^-l)j 



(5 + 4j(r?- 1) 



2 ' 



J4).n^_o , l + 3r] + j(^' 



g^'\Q) = 2 + —^^±^-^^^^. (11) 
(l+j(^-l) 



An interesting result here is that the above correlation functions behave dif- 
ferently for even (j = 1/4) or odd (j = 3/4) photon numbers. For instance, 
c;(2)(o) = (3 + r])/2 and c/(^)(0) = 2 + 4(3 + 12r/ + rf)/{?> + rff when j = 1/4, 
while ^(2)(o) = 2(1 + 3r/)/(l+ 7^)2 and ^(^)(0) = 2 + 4(3 + 12r/ + 3?72)/(l + Sr/)^ 
when j = 3/4. Particularly, for even or odd photon number states, the second- 
order coherence function g^'^\Q) will be linearly or inversely proportional to 
f] when r] increases. Depending on the steady-state behaviors of the photon 
correlation functions, one can distinguish between the nonlinear oscillator's 
odd and even photon number states. Thus, the photon coherence functions 
are a convenient tool to determine the parity of the photon number of the 
state |j, m) which corresponds to su(l,l) algebra. 

In Fig. ([2]) we show the dependence of the second-order coherence function 
g'^'^\Q) versus the detuning Su and different values of j. These behaviors are 
explained as follows: for an even number of photons describing the state | j, m) , 
i.e. j = 1/4, the mean-photon number (n) will be below unity (for n = 2) and 
the normalized second-order coherence function increases accordingly, show- 
ing super- Poissonian photon statistics {g^'^\0) > 1). Conversely, for an odd 
number of photons, i.e. j = 3/4, the mean-photon number (n) will be above 
unity (or near unity) and the second-order coherence function decreases, re- 
vealing near Poissonian (^''•^■'(O) ~ 1) or even sub- Poissonian (^''■^^(O) < 1) 
photon statistics and a single-photon state can be created here. Note that the 
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fourth-order coherence function g^^\0) approximately behaves as g^'^\0), but 
with a different magnitude. 



4 Summary 

In summary, we described a scheme capable of cooling an oscillator coupled 
to an externally pumped artificial atom (a Josephson flux qubit) and damped 
nonlinearly through interaction with its environmental thermal reservoir. Un- 
der certain conditions, the oscillator and the qubit exchange two-photons, 
allowing us to describe their quantum dynamics using the su(l,l) algebra. If 
the qubit 's dynamics is faster than that of the nonlinear oscillator, the cool- 
ing of the oscillator's degrees of freedom occurs when controlling the qubit 
quantum dynamics. Evaluating the second-order photon correlation function 
(or higher-order correlations), one can distinguish between even and odd pho- 
ton number states characterizing the oscillator. By adjusting the parameters 
involved, one can create a nonclassical field state with sub-Poissonian photon 
statistics. This will allow us to obtain a single-photon state of the nonlinear 
oscillator. 



A Appendix 

In this Appendix we obtain the equation ([6]). We start by indicating the Hamil- 
tonian Hf describing the interaction of the environmental bath with the non- 
linear oscillator, i.e. the Eq. (jlj): 

Hf = hub^b + 2nx{b^(3- + (3+b). 

In the Born-Markov approximation one can eliminate the bath operators. For 
doing this we define an operator Q f which belongs to the oscillator's subsystem 
and satisfy the following equation of motion: 

j^iQj) = mmp-,Qf]) + {[Qf,f3^]b)}. (A.i) 

The formal solution of the Heisenberg equation for b^ is: 



9 



v — 2uc + ix 
V — 2uJc + IX 



(A.2) 



with h{t) = [6^(t)]+. Substituting Eq. ( ICTl in Eq. ( ICT) and using the Bogol- 
ubov lemma l26l 



where U is an arbitrary operator belonging to the oscillator subsystem together 
with the identity Tr{f^Qf{t)pf{0)} = Tr{|p/(t)Q/(0)} one arrives at Eq. 

There k = , , and x = 
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